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We investigate a selection-mutation model for the dynamics of technological innovation, a spe- 
cial case of reaction-diffusion equations. Although mutations are assumed to increase the variety 
of technologies, not their average success ("fitness"), they are an essential prerequisite for innova- 
tion. Together with a selection of above-average technologies due to imitation behavior, they are 
the "driving force" for the continuous increase in fitness. We will give analytical solutions for the 
probability distribution of technologies for special cases and in the limit of large times. 

The selection dynamics is modelled by a "proportional imitation" of better technologies. However, 
the assessment of a technology's fitness may be imperfect and, therefore, vary stochastically. We will 
derive conditions, under which wrong assessment of fitness can accelerate the innovation dynamics, 
as it has been found in some surprising numerical investigations. 

PACS numbers: 89.65. Gh,47.70.-n,89.20.-a,47.54.+r 



Introduction. Innovation dynamics has not only been 
studied by economists and social scientists (see Refs. in 
0), but also recently attracted an increasing interest by 
physicists [J 0,0. Ebeling et al, for example, investigate 
the competition of an innovation with already established 
technologies 0,IH- 

In a wider perspective, the selection-mutuation equa- 
tions discussed in this manuscript are a special case of 
reaction-diffusion equations, where technological imita- 
tion processes are analogous to chemical reactions and 
diffusion originates from mutation processes. Due to their 
non-linearity, reaction-mutation equations are known to 
display interesting pattern-formation phenomena such as 
(spiral) waves |y,ul or Turing patterns pjQ- However, 
this subject will not be the focus of this paper. Instead, 
we will try to understand some surprising observations 
in a microsimulation model of innovation dynamics by 
Saam 

Our innovation model assumes that the success of a 
company's technology can be expressed by a "fitness" 
value x. The technologies of superior companies are im- 
itated (copied) with a certain probability per unit time 
specified later. Moreover, all companies perform innova- 
tive activities, which are modeled by random mutations 
of the currently applied technologies. Assuming a Brown- 
ian motion for the mutation of fitness values, they change 
from some value x at time t to some other value z at time 
t + At with a Gaussian transition probability 

p -[z(t+At)-x(t)] 2 / (28V At) 

p(z(t + At)\x(tj) = ===== . (I) 

VzttOv At 

Here, v is the frequency, v At the number and Vd the 
standard deviation (average size) of mutations. Accord- 
ing to Eq. Q), mutations may be beneficial, but they 
may also decrease the fitness of technologies. Therefore, 
the observed increase in the average fitness would have 



to be a consequence of the assumed imitation of superior 
technologies. 

Patents on the leading technologies could suppress 
such imitation and, therefore, potentially reduce the 
speed of innovation, i.e. the increase in the average fit- 
ness of technologies. However, under certain conditions, a 
misperception of fitness values Q] seems to neutralize or 
even over-compensate such a deceleration effect 0. Our 
paper will try to achieve an analytical understanding of 
these interesting findings. 

Formulation of the Innovation Model. For the sake of 
analytical treatment, let us assume a continuous spec- 
trum of possible fitness values x. Moreover, let P(x, t) dx 
represent the probability with which the applied tech- 
nology of a company has a fitness value between x and 
x + dx at time t. In addition, we will assume that compa- 
nies imitate technologies with a higher fitness z > x and 
that the imitation rate increases proportionally with the 
related increase z — x in success. If z — x is negative, the 
imitation rate should be zero. Then, given the presently 
applied technology has the fitness x, the respective imita- 
tion rate w± {z\x,t) to copy another technology of fitness 

z, is dJGjl 



wi(z\x, t) = XP{z, t) max(z — x, 0) 



(2) 



as the imitation rate is proportional to the occurence 
probability density P(z,t) of the imitated technology z. 
The parameter A > has the meaning of an imitation 
frequency. 

In order to calculate the change of the occurence prob- 
ability density P(x, t) due to imitation processes, we have 
to insert w(z\x, t) = wi(z\x, t) into the continuous master 
equation 0,0, El 



dP(x, t) 
dt 



= / dz [w(x\z,t)P(z,t) — w(z\x,t)P(x,t)) , 



(3) 



2 



which leads to 
dP(x, t) 

dt 



A / dzP(x,t)P(z,t) 

x [max(x — z, 0) — max(z — x, 0)] (4) 



= XP(x,t) 



x — dz zP(z, t) 



(5) 



av. fitness (x)t 



because of the normalization condition J dzP(z,t) = 1. 
Equation is related to the selection equation from 
evolutionary biology [Hi. Il4| and to the game-dynamical 
equations |10llllUl5| . In order to take into account new 
inventions, with 



w(z\x. t) = lim 
At^o 



p(z(t + At)\x(t)) 
At 



■w\{z\x,t) (6) 



we will additionally assume random transitions to other 
technologies z with the rate Wo(z\x,t) = w(z\x,t) — 
w\{z\x,t). According to this, we simply treat inventions 
as unbiased random mutations, which may either increase 
or decrease the fitness of the pursued technology, as in 
complex (technological) systems it is sometimes hard to 
judge before, whether some change will cause an improve- 
ment or not. Although the assumption of unbiasedness 
may be questioned, it is favorable to gain non-trivial con- 
clusions. The resulting equation for the temporal change 
of P(x, t) can be approximated by the following Fokker- 
Planck type of equation [16| : 



dP(x,t) 
dt 



XP{x,t) 



x — I dz zP{z, t) 



D 



d 2 P(x,t) 
dx 2 



(7) 

The quantity D = v9/2 has the meaning of a diffusion 
coefficient. For the mean fitness (x) t one can derive the 
differential equation flTL f]~3| 



d(x) 
dt 



AVar t (x) 



(8) 



while for the variance Var t (ir) = ((x — (x)t) 2 )t, we find 



dVa,r t (x) 
dt 



= \({x-(x) t ) 3 ) t + 2D. 



(9) 



The term ((x — (x) t ) 3 )t disappears if the distribution 
P(x, t) has a vanishing skewness 7, e.g. for a Gaus- 
sian distribution. Equation JHJl shows that the diver- 
sity Var t (x) is a key factor for a high innovation speed 
d{x)t/dt. 

Equation J7J) can be viewed as a special case of 
reaction-diffusion equations 0, 0, where the diffusion 
term originates from mutations and the non-linear re- 
action term from imitative pair interactions Eaini. it 

is known to have a complicated formal solution, which 



is uniquely determined by the initial condition [TTL Il3| . 
However, this does not help us, here. We are rather look- 
ing for a qualitative understanding of the dynamic solu- 
tion and, where possible, for explicit analytical expres- 
sions. In the limiting case v = = D of no muta- 
tion and no diffusion, all companies will eventually im- 
itate the technology with the highest fitness x max , i.e. 
limt_>oo P{x, t) = 6(x — x max ), where S(x) denotes Dirac's 
delta function. In the limiting case A = of no imitation, 
Eq. J7J) is just a linear diffusion equation, and its solu- 
tion is a superposition of Gaussian distributions with a 
linearly- in-time increasing variance 2Dt: 



P(x,t) 



d z P ( Z ^ p -(x-zf/(±Dt) 



(10) 



If we initially have a two-point distribution P{x, 0) = 
P\5{x — Xi) + P2^(x — X2) with px + P2 = 1, both peaks 
will become broader due to diffusion similarly to Eq. I jlOp . 
At the same time, however, imitation tends to favour the 
peak around the superior technology's fitness, while the 
peak around the inferior technology's fitness will eventu- 
ally disappear. After some time, one finds a unimodal, 
almost Gaussian distribution, if A and D are constant, 
see Fig. ^ 
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FIG. 1: (Color Online) (a) Computer simulation of Eq. J7J 
with D = 1 and A = 0.1. Starting with a two-peak initial dis- 
tribution P(x, 0), one can see that the superposition law does 
not apply due to the non-linearity of our innovation model. 
Nevertheless, after some time (t > 5) we find the Gaussian 
distribution Jilt , (b) After the initially developing skewness 
7(t) = ((a;-(a;} t ) 3 } t /Var t (a;) 3/2 has vanished (t > 5), the vari- 
ance Var^x) grows linearly and the mean fitness {x) t quadrat- 
ically, in accordance with Eqs. il Hil and l)17ft . 



Therefore, let us investigate how a Gaussian distribu- 
tion develops in time. We will show that the non-linear 
Eq. 10 has, then, a special solution of the form 



P(x,t) 



[x-f(t)] 2 /[2g(t)] 



(11) 



where f(t) denotes the mean value and g(t) the variance 
of the fitness values x. Applying the product, quotient, 



3 



and chain rules of differential calculus, we find 
dP(x,t) _ 2[x-f(t)} 



dx 
d 2 P(x,t) 
dx 2 



dP(x,t) 

m 



<-P(x,t) 



(12) 



= >" + ir*" (13) 



1 ^£ 



-l*P(z,i) + 



2 3(t) 2 



P(i,t). 



(14) 



Comparing the expression for dP(x,t)/dt with the one 
for d 2 P(x,t)/dx 2 shows that Eq. ifTTl) implies 



9x 2 



<?(*) 



Considering ((x - (x) t ) 3 )t = 0, f(t) = (x) t , g(t) = 
Var t (x), and Eqs. JBJ), ©, this partial differential equa- 
tion indeed agrees with Eq. (0. Therefore, Eq. 1)11(1 is 
the unique solution if the initial condition is a Gaussian 
distribution. The existence of such a generalized diffusion 
solution is quite surprising in view of the non-linearity of 
Eq. J7J. This reminds one of the exact solution of the 
non-linear Burgers equation |17|. 

In the previously assumed case D = const, of constant 
mutation activity, we obtain 



g(t) = g(0) + 2Dt, 

f(t) = f(Q) + Xg(0)t + XDt 2 



(16) 
(17) 



However, for D(t) = Co<?(i), i.e. if the overall mutation 
activity would be proportional to the variety of existing 
technologies, we would find an exponential growth 

g(t) = g(0)e 2C ^, (18) 
/(*) = (/(°)-^)+^ 2C0t - ( 19 ) 

Finally, for D(t) = Ci/(t), i.e. if the overall mutation 
activity were proportional to the average fitness (and the 
related profits), we would expect 

VXg(0) 



f(t) = f(0)cosMV2XCit) 



sinh( v / 2ACii) . 



In summary, starting with a normal distribution of 
technologies x, the solution of Eq. Q stays normally dis- 
tributed, but the average f(t) and variance g(t) depend 
significantly on the diffusion coefficient D and, hence, on 
the mutation activity, i.e. v and 9. If the diffusion coeffi- 
cient stays constant, the average fitness grows quadrati- 
cally in time, while it tends to grow exponentially, if the 
diffusion coefficient increases proportionally to the mean 



value f(t) or the variance git) of the fitness of technolo- 
gies. An exponential growth is probably the more realistic 
scenario (see, for example, Moore's law). 

Imperfect Evaluation of Fitness. Let us now discuss 
a generalization of the above model, assuming that the 
fitness is not assessed exactly. For this, let £ denote the 
difference between the perceived fitness and the actual 
fitness a; of a technology, and z + £ the perceived fitness 
of a technology with fitness z. In this case, we have the 
formula 



w\{z\x,t) — A max(z + £ — x — £ , 0)P(z, t) 



(21) 



and study the probability density P'(x,£,t) = 
Px(QP(%,t), wnere Px(0 denotes the conditional prob- 
ability P((,\x) that the error in the fitness estimate is £, 
given the actual fitness is x. The following normalization 
relation applies: J d£P x (£) = 1- Therefore, the resulting 
differential equation for the distribution P(x, t) of fitness 
values becomes 



(15) dP(x,t) 



dt 



A / dz I d4 I d(P x (Z)P(x,t)P z (QP(z,t) 

d 2 P(x,t) 



x(x + £- z-Q+D- 



dx 2 



J d?P x (Z). (22) 



The solution of this equation depends crucially on the 
distribution P x {^) and can, in general, not be analyti- 
cally determined. We will, therefore, restrict to deriving 
a generalization of Eq. (JHJ in order to see, how imperfect 
evaluation of a technology's fitness affects the temporal 
evolution of the average fitness, if at all. Applying the 
method of partial integration to the diffusion term twice 
and defining the mean value of a function G{x, £, z, £, t) 
by 



(G) t 

we obtain 

d(x) t 
dt 



dx J d£ J dz J dC, G(x,£, z,C, i) 
xP(x,t)P x (£)P(z,t)P,(Q, (23) 



dx x 



dP{x,t) 
dt 



(24) 



X I dX I ^ I dZ I d ^ P ^ t ) Px< ^) P ^ z ^ 

(25) 



xP z (()(x 2 +xH — xz — x() 
d 2 P(x,t) 



D / dx 



dx 2 



= A [(x 2 ) t + (xOt - (xz) t - (x() t ] 
+ D J dx^P(x,t). 



(26) 



=o 



4 



With (z) t = {x) t , (Ot = {Qt, {xz) t = (z)t<2>t, 
Cov t (x,C) = 0, Cov t (x,f) = (xQ t - (x) t (Ot, and analo- 
gous relationships, we find 

^ = A [Var t (x) + Cov t (x,£)] . (27) 

Compared to Eq. imperfect perception of fitness 
results in the additional term Covt(cc,£). Accordingly, 
misperccption does not have any influence on the dynam- 
ics, if the deviations £ are statistically independent of x 
or P x (£) is not a function of £. If Cov t (x, £) < 0, i.e. if 
leading technologies are systematically underestimated, 
while inferior technologies are overestimated, mispercep- 
tion slows down the increase of the average fitness. How- 
ever, misperception will speed up the evolution of better 
technologies, if Cov t (a;,£) > 0, i.e. if leading technolo- 
gies are overestimated. This is actually the case in the 
simulation model by Saam [9(, as it assumes that only 
the technology of the firm with the best perceived fit- 
ness is copied, while other firms are never imitated, even 
if their technologies are better than the own one. The 
same simulation study has shown that an overestimation 
of the apparently best technology can even compensate 
for inhibitory effects of patents, which were assumed to 
suppress the imitation of the actually fittest technology 

Summary. In this paper, we have analytically inves- 
tigated the dynamics of a selection-mutation model for 
technological innovations, which can be viewed as a spe- 
cial reaction-diffusion model, but without emergent pat- 
tern formation based on a Turing instability 0, Q . In 
our case, the imitation of better technologies plays a role 



analogous to chemical reactions. Imitation has an effect 
comparable to the selection of above-average strategies in 
similar equations of biological evolution. Inventions have 
been modeled as unbiased, random mutations. They may 
improve the fitness of technologies or deteriorate them. 
As a consequence, the observed increase of the average 
fitness requires a selection process. However, a selection 
process alone also does not cause a persistent increase in 
the average fitness: If the mutation frequency v and the 
diffusion coefficient D are 0, the fitness of all technologies 
converge to the highest initial fitness and stay there. Only 
if both, random mutations of technologies and imitation 
are combined, we have a steady growth of the average fit- 
ness, so that after some time even a bad technology will 
be replaced by another one which is much better than the 
initially leading technology. Interestingly, the speed of the 
average increase in fitness is proportional to the variance 
in the fitness of technologies. Therefore, copying other 
companies' technologies alone does not support a persis- 
tent innovation trend. Instead, diversity is the "motor" 
or "driving force" of innovation. 

Finally, we have investigated the effect of mispercep- 
tion of the fitness of technologies. It turned out that mis- 
perception can be neutral, but it can also speed up or slow 
down technological evolution. Our results indicated that 
the average fitness will grow faster, if the leading tech- 
nologies are systematically overestimated and the fitness 
of inferior technologies is underestimated. Therefore, ex- 
citement for new technologies can speed up innovations 
even without higher investments, just because of a bias 
in the perception of fitness, i.e. a bias in the imitation 
behavior of superior technologies. 
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